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ABSTRACT 

We present a rest-frame ultraviolet morphological analysis of 108 z ~ 2.1 Lyman Alpha Emitters 
(LAEs) in the Extended Chandra Deep Field South (ECDF-S) and compare it to a similar sample of 
171 LAEs at z ~ 3.1. Using Hubble Space Telescope (HST) images from the Galaxy Evolution from 
Morphology and SEDs survey, Great Observatories Origins Deep Survey, and Hubble Ultradeep Field, 
we measure size and photometric component distributions, where photometric components are defined 
as distinct clumps of UV-continuum emission. At both redshifts, the majority of LAEs have observed 
half-light radii < 2 kpc, but the median half-light radius rises from 1.0 kpc at z = 3.1 to 1.4 kpc at 
z = 2.1. A similar evolution is seen in the sizes of individual rest-UV components, but there is no 
evidence for evolution in the number of multi-component systems. In the z = 2.1 sample, we see clear 
correlations between the size of an LAE and other physical properties derived from its SED. LAEs are 
found to be larger for galaxies with higher stellar mass, star formation rate, and dust obscuration, but 
there is no evidence for a trend between equivalent width and half-light radius at either redshift. The 
presence of these correlations suggests that a wide range of objects are being selected by LAE surveys 
at z 2, including a significant fraction of objects for which a massive and moderately extended 
population of old stars underlies the young starburst giving rise to the Lya emission. 
Subject headings: cosmology: observations — galaxies: formation - galaxies: high-redshift - galaxies: 
structure 



1. INTRODUCTION 

The majority of galaxies in the local universe 
lie on the Hubble sequence (IHubbld [T936h . a con- 
tinuum that runs between red, passively-evolving 
with a compact spheroidal component and 
star-forming disks with exponential light 
This morphological sequence is seen 



galaxies 
gas-rich, 
profiles. 

clearl y out to intermediate redshift s (z 



e.g., Glazebrook et al.lll995l: Ivan den Berah ct al 
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Griffiths et al.lll996l: iBri nchmann eTal. 1998; Lillv et atl 



19981: iSimard et al l 11999; van Dokkum ct al. I200C- 
Stanford et al.ll200^ iRavindranath et al.ii2004i ). beyond 
which the majority of g alaxies appe ar clumpy and irreg- 
ular (e.g., 'Giavalisco ct al. 1996; Lo wenthal et al.l 119971: 
Dickin son' 2000: van den Bergh .20011: iPapoyich et al.l 



2005; Conschc e'et al.ll2005l : iPirzkal et al.il2007| ) and are 

difficult to place into existing classification schemes. At 

^ Centro de Astro-Ingenieria, Pontificia Universidad Catolica 
de Chile, Santiago, Chile 



2 ^ 2 ^ 3.5, these galaxies, the majority of which are 
actively star-forming, have sizes ranging from < 1 kpc 
to ~ 5 kpc, with the largest ofte n exhibiting multi- 
ple ph otometric components (e.g., Bouw ens et al.|[2004l : 
iRavindr anath et al. 2006; Oesch ct al. 2009D. 

One of the best-studied classes of galaxies at high red- 
shift is the Lyman Alpha Emitter (LAE). At z ~ 3 - 4, 
LAEs are widely believed to be actively star-forming, 
low in both stellar and dark m atter mass, and rela- 
tive ly dust-fre e (e.g. , lCowie fc Hu 1998; Vcncn ians et al.l 
[200,4 iGawiser"et~aLll2007n . The continuum morpholog- 
ical properties of LAEs vary from object to object, but 
the majorit y are compact (C > 2.5) and have size s 
< 1.5 kpc (iVenemans e t aU '2005"; 'Pirzkal et al.l 12007 



Overzier et al.l 120081: iTan iguchi et al. 2009; Bon d et al.l 
20091 : iGronwall et al.ll2010l) . Multi-component or clumpy 
LAEs make up ~ 20 — 45% of the population and 
typically have morphologies that are qualitatively sim- 
ilar to other star-forming galaxies at high redshift. 
The emission-line morphologies of LAEs are difficult 



2 



Bond et al. 



TABLE 1 
LAE Samples 



Name Size HST Coverage^' GEMS GOODS HUDF Reference 



z2Guaita 250 193 175 32 6 1 

z'2EWComplete'' 130 108 97 19 6 1 

z3Gronwall 154 116 97 29 3 2 

z3Ciardullo 62" 55" 47" 12" 0" 3 



References. — (1) Guaita et al. 2010; (2) Gronwall et al. 2007; (3) CiarduUo et 
al. 2011 

^ Counts the total number of objects covered by at least one of GEMS, GOODS, or 
HUDF 

z2EWcomplete is a subsample of z2Guaita 
" Counts only those objects that are not already counted in z3Gronwall 



to study becau s e of t he l ong exposure times requ ired, 
but iBond et all (poToh and iFinkelstein eTall (poTol ) find 
them to be spatially compact (re < 1.5 kpc) at z — 3.1 
and z — 4.4, respectively. There may be evidence of an 
additional extended and diffuse component of Lya emis- 
sion in LAEs at z = 2.3 (iNilsson et al.ll2009l) and z = 4.4 
(IFinkelstein et all [20Tol ). but iBond etalTpoTol ) found 
no evidence for this component in LAEs at z = 3.1, so 
the question r emain s open. Using ground-based imaging, 
iSteidel et all (|201lD found extended Lya emission out to 
~ 80 kpc around a stack of LAEs identified by the Ly- 
man break technique, but it is unclear to what extent 
these highly extended halos contribute to the emission- 
line morphologies of typical LAEs on kiloparsec scales. 

The Mul tiwavelength Survey by Yale-Chile (MUSYC, 
lGawiser~et al. 2006) has obtained multiwavelength imag- 
ing and spectroscopy of 1.2 degree^ of sky in four 
fields, including the Extended Chan dra Deep Field - South 
(ECDF-S). As part of this survey, iGuaita et all (|2010f ) 
(hereafter, GulO) used broadband and 3727 A narrow- 
band imaging of the ECDF-S to identify a large, unbiased 
sample of LAEs at z = 2.1. The authors found the LAEs 
in this sample to be weakly clustered, with a bias fac- 
tor, b ^ 1.8, that is similar to that expected from the 
progenitors of prcscnt-day L* galaxies. An analysis of 
the b roadband o ptical and infrared colors of a "stacked" 
LAE (jCuaita et ai.ll201]] ) suggests that their median stel- 
lar ma sses (~ 4 x 10^ Mfr^) are similar to those of z = 3.1 
LAEs (jGawiser et al.ll2007l ). although there is a subset of 
IRAC-bright objects that is ~ 10 times more massive, on 
average, and exhibits no n-negligible amounts of dust red- 
denin g, £'(B-V) - 0.4 (|Lai et al.ll200llAcquaviva et al.l 

[mil). 

In this paper, we study the rest-UV continuum 
morphologies of the GulO sample of LAEs and com- 
pare them to those seen in LAEs at z = 3.1 
(jBond eta L 2009; Gronwal l et al.ll2010( ) using the sam- 
ples of iG^onwa ll et al. (20o3' and a new sample from 
ICiardullo et aH (12011 7 We use images taken by the Ad- 
vanced Camera for Surveys (ACS) and obtained as part 
of the Galaxy Evolution from Morphology and SEDs sur- 
vey fCEMS.iRbc et al.l l2004f). Great Observatorie s Ori- 
gins Deep Survey ('GOODS. ICiavalisco et aIll2004D. and 
Hubb le Ultradeep Field survev (HUDF. te eckwith et al.l 
120061) . In what follows, we fit each photometric compo- 
nent separately and give quantitative size measures for 
both the individual components and the LAE system as 
a whole. 



In § [2] and [3l we summarize the data and describe the 
analysis techniques used in our comparative study. In 
§ 131 we give the photometric properties of each z = 2.1 
LAE system and its components and compare our re- 
sults to those found for z = 3.1 LAEs. Finally, in § |5l 
we discuss the implications of our findings for the phys- 
ical nature of LAEs as a function of redshift. Through- 
out this paper, we will assume a concordance cosmology 
with Hp = 71 km s" ^ Mpc'S = 0.27, and fl^ = 0.73 
(jSpergel et al.ll2007[ ). With these values, 1" = 7.75 phys- 
ical kpc at z = 3.1 and 1" = 8.43 physical kpc at z = 2.1. 

2. DATA 

We use three LAE samples in our analysis, includ- 
ing the equivalent- width-complete subsample of z — 2.1 
LAEs identi fied by GulO, the flux -limited sample of z ~ 
3.1 LAEs of iGronwall et al.l (|2007l z3Gronwall), a nd the 
flux-li mited sample of z ~ 3.1 LAEs of fCiarduUo et all 
IpOnl z3Ciardullo), all in the ECDF-S. The Veoe-band 
cutouts and morphol ogical prope r ties o f z3Gronwall are 
already published in IBond et al.l (|2009L hereafter B09). 
All samples are summarized in Table [T] and described 
below. 

The initial sample of 250 z = 2.1 LAEs (hereafter, 
z2Guaita) was selected to have a monochromatic flux 
of -FLyQ > 2 X 10"^'' ergs cm~^ s~^, and rest-frame 
Lya equivalent width of EW > 20 A. For some anal- 
yses, the authors made a further cut on Lya luminos- 
ity, Llyq > 1-3 X 10"^^ erg s~^, in order to remove a 
bias in their sample against high-EW objects. We use 
this "equivalent-width complete" subsample (hereafter, 
z2EWcomplete) of 130 z = 2.1 LAEs for the majority of 
our morphological analyses. Excluding objects within 40 
pixels of the edge of an image and cutouts with clear im- 
age defects, there are 108 z = 2.1 LAEs in z2EWcomplete 
that are covered by HST broadband imaging surveys. 

The two z ~ 3.1 LAE samples used here were both 
selected to have monochromatic fluxes, FLya > 1-5 x 
10~^^ ergs cm~^ s~^, and rest-frame Lya equivalent 
widths, EW > 20 A. The z3Ciardullo sample covers 
the same field as z3Gronwall, but identifies LAEs using 
a narrow-band filter shifted 20 A to the red relative 
to the z3Gronwall filter. There remains some overlap be- 
tween the two filters, however, so ~ 50% of the objects 
identified in z3Ciardullo also appear in z3Gronwall. We 
find 55 additional LAEs in z3Ciardullo that have HST 
coverage, complementing the 116 z ~ 3.1 LAEs stud- 
ied in B09. Note that we also exclude objects in these 
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samples that are wit hin 1" of any X-ray s ource in the 
expanded c atalogs of iLehmer et alj ()2005D , iVirani et alj 
((200l) . and lLuo et all (I2008D . 



2.1. GEMS 

Of the surveys used in our study, the GEMS survey 
covers the largest area, consisting of ^ 800 arcmin^ of 
the ECDF-S and 63 ACS pointings in the F606W (Veoe) 
and F814W (zg^) hhers. The depth of GEMS is rela- 
tively uniform across the field, with a Sct detection limit 
of rriAB = 28.3 for Vsnfi- band point sources. S urvey 
images were multidrizzled (jKoekemoer et all 120021 ) to a 
pixel scale of 30 mas. The cutouts for the LAEs in the 
z = 2.1 and z = 3.1 samples are shown in Figures [T] and 
[2] (see end of article) , respectively. 

2.2. GOODS 

The southern half of the GOODS survey is a subregion 
of the Extended Chandra Deep Field-South, observing ^ 
160 arcmin^ of sky. The GOODS observations included 
HST/ ACS imaging in the F435W (B435), Vgoe, F775W 
(I775), and Z850 filters. Although the Vgoe imaging depth 
varies across the GOODS area, a typical 5cr detection 
limit for point sources is toab = 28.8. As in GEMS, 
all images were multidrizzled to a pixel scale of 30 mas. 
The cutouts for the LAEs in the 2; = 2.1 and z = 3.1 
samples are shown in Figures |3] and 2] (see end of article), 
respectively. 

2.3. HUDF 

The Hubble Uhra-Deep Field (HUDF) reaches a Vgoe- 
band 5a point source depth of toab = 30.5 and includes 
HST/ ACS imaging in the B435, Vgoe, I775, and zgso fil- 
ters. The survey covers 11 arcmin^ of sky and has a 
multidrizzled pixel scale of 30 mas. The cutouts for the 
LAEs in z2EWconiplete are shown in Figure [5] (see end 
of article). 

3. METHODOLOGY 

3.1. Visual classification 

We have classified all of the 193 objects in z2Guaita 
with HST coverage according to their morphology in the 
Vgoe cutouts. The morphological classifications are de- 
termined by eye, taking into account an object's size 
and appearance, where the classes are chosen to sepa- 
rate objects with unusual physical conditions (e.g., ap- 
parent major mergers), as well as possible interlopers or 
contaminants. The four classes are given below. 

• LCo - Compact objects. Much like the majority of 
LAEs at z = 3.1, these objects are compact with 
little or no evidence for extended substructure. 

• CEx - Clumpy-extended objects. They appear sim- 
ilar to the most extended high-redshift star-f orming 
galax ies (e.g., clump-clusters, rElmcgrcc n et al.l 
l200i ). with emission appearing in disconnected 
clumps. 

• DEx - Diffuse-extended objects. These objects are 
large enough to be at least marginally resolved from 
the ground, with emission dominated by a diffuse 
component. If they really are at z = 2.1, they are 
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Fig. 6. — Example cutouts of objects that are visually classified as 
diffuse-extended (DEx, bottom row), clumpy-extended (CEx, mid- 
dle row), and compact (LCo, top row). Galaxies in the DEx class 
are excluded from some of our analyses because of the high prob- 
ability that they are coincident with low-rcdshift contaminants. 
These cutouts are 2" 4 on a side. 
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Fig. 7. — Distribution of z = 2.1 LAEs in rest-frame equiva- 
lent width vs. Lya luminosity, with symbols indicating the visual 
morphological classification of the candidate LAE. Objects with 
classifications of DEx, CEx, LCo, and NoD are plotted as stars, 
filled squares, filled triangles, and open circles, respectively. The 
horizontal solid line is the 20 A rest-frame equivalent width cut- 
off of the GulO survey and the dotted vertical line indicates the 
log(LLya) > 42.1 luminosity cut used to make the z2EWcomplete 
sample. The dashed curve delimits the selection region based on 
the 5o- limiting narrow-band magnitude of the GulO survey. 

very massive systems, but many are likely to be 
contaminants or low-redshift interlopers coincident 
with a high-redshift LAE. 

• NoD - Non-detected objects. There is no discern- 
able source in the Veoe cutout. 

Examples of the LCo, CEx, and DEx classes are shown in 
FigurelH Of the objects in z2Guaita with HST coverage, 
a strong majority (71%) are classified as LCo, with the 
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Fig. 8. — Distribution of SExtractor detections in the Veoe-band 
GEMS cutouts as a function of distance from the ground-based Lyo 
centroid. We classify as LAE components any detections within the 
0''65 (22-pixeI) selection radius, drawn in black. 

CEx, DEx, and NoD classes making up 10%, 16%, and 
3% of the sample, respectively. 

The candidate LAEs are plotted in equivalent width vs. 
Lya luminosity in Figure [71 where symbol type indicates 
the visual classification. GulO selected LAEs accord- 
ing to their narrow-band magnitude rather than their 
Lya luminosity, leading to a bias against high-equivalent- 
width LAEs with log(Lya) < 42.1. Therefore, the ma- 
jority of our analyses are performed on the subsample 
with log(Lya) > 42.1, also referred to as z2EWcomplete. 
The use of this subsample is further motivated by the 
small fraction (6%) of objects in the DEx class with 
log(LyQ;) > 42.1, as these objects may be dominated 
by low-redshift contaminants or interlopers within the 
selection radius. 

3.2. Source extraction and aperture photometry 

In B09, we describe in detail a data analysis pipeline 
optimized for measuring the photometric properties of 
"clumpy" or irregular star-forming galaxies. We will pro- 
vide only a brief summary here. 

For the zSCiarduUo sample, we followed B09 and ex- 
tracted 80 X 80 pixel (2'.'4 x 2'.'4) cutouts from the GEMS, 
GOODS, and HUDF images at the ground-based po- 
sition of each LAE in our sample. Our final sam- 
ple includes only those LAEs with full survey cover- 
age in the cutout region. We identified all sources in 
each cutout using SExtractor (Bcrtin & Arnouts 1996) 
with extraction parameters DETECT_MINAREA= 30 and 
DEBLEND_MINCONT= 0.06, fitting and subtracting a uni- 
form sky from each cutout. In order to determine the 
rest-UV centroid of each LAE system, we again run SEx- 
tractor on each cutout, now with DETECT_MINAREA= 5, 
and take the flux-weighted mean position of the result- 
ing detections. We follow the same procedure for the 
z2EWcomplete sample, with the exception that we ex- 
tract somewhat larger cutouts, 100 x 100 pixel (3" x 3"), 
to provide full coverage for some of the more extended 
objects in the z2Guaita sample. 

For the z3Ciardullo sample, we follow B09 and se- 



lect sources within rgei = 0'.'6. For z2EWcomplete, in 
order to account for the difference in the angular di- 
ameter distance between z = 3.1 and z = 2.1, we set 
Tsei = 0'.'65. In Figure m we plot the distribution of SEx- 
tractor Veoe-band detections as a function of angular dis- 
tance from the ground-based Lya centroid for the objects 
in z2EWcomplete. Based on the density of sources out- 
side the selection region, we estimate that at most ^ 12% 
and ~ 14% of all LAEs will have an interloper within rgoi 
in z3Ciardullo and z2EWcomplete, respectively. The ac- 
tual contamination rate will likely be less than this be- 
cause the presence of an interloping source within ^ 1" 
will decrease the apparent equivalent width of the LAE 
in ground-based imaging and make it less likely to exceed 
the EW cutoff for the surveys. 

Following B09, we define the photometric centroid to 
be the flux-weighted mean position of the detections 
within Tsoi. Aperture photometry and half-light radii are 
computed within fixed apertures with radius, Tap = rsc\, 
using the IRAF routine PHOT and centered on the cen- 
troid determined from the SExtractor runs. 

3.3. Monte Carlo simulations 

The HST/ ACS images used in this study contain 
strong pixcl-to-pixe l correlat ions as a result of the driz- 
zling process (Kock emoer et al. 2002), so we estimate the 
uncertainties on the continuum magnitude and half-light 
radii using Monte Carlo simulations. We perform a to- 
tal of 10^ simulations for each survey, each time placing 
a Gaussian profile at a random position on the Veoe- 
band image. The simulated sources have Gaussian pro- 
files with a dispersion between 2 and 12 pixels (between 
0'.'06 and 0'.'36) and their photometric properties are com- 
puted within a fixed O'.'G aperture. We define the scatter 
on a photometric measurement in terms of median statis- 
tics: 

< = 0.748[Q3(a;)-Ql(x)], (1) 

where Ql and Q3 are the first and third quartiles, re- 
spectively. In the limit of a well-sampled Gaussian dis- 
tribution, this quantity approaches the square root of the 
variance. 

The uncertainty on the fixed-aperture magnitudes can 
be estimated using the weight images provided with each 
survey, where the pixel-to-pixel flux uncertainty is given 

by fj/ = y'lJW, with the mean being taken over pixel 
weight values, W, in each 0'.'6 aperture. When we ana- 
lyze the scatter in our simulations, we find a linear rela- 
tionship between the fractional uncertainty in half-light 
radius and the fractional uncertainty in flux within the 
aperture (see Figure ^ . For the Gaussian model profiles 
used in our simulations, a least-squares fit yields, 

^ = 0.54^. (2) 

re /606 

As shown in the figure, this relationship holds for all of 
the surveys used here (GEMS, GOODS, and HUDF). 
Furthermore, Figure [TOl demonstrates that the fractional 
uncertainty in half-light radius is independent of half- 
light radius over the range typical for LAEs 0'.'06 — 
0'.'15). 

As noted in B09, there is a systematic tendency to over- 
estimate the fixed-aperture half-light radius at very faint 
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TABLE 2 
2.1 LAE Properties 



Fig. 9. — Fractional uncertainties in the half-light radius as a 
function of the fractional uncertainty in flux through the Vgoe fil- 
ter for a scries of simulated Gaussian-shaped light profiles placed 
at random positions on the survey images. Flux uncertainties are 
determined from the survey weight images and symbol type indi- 
cates the survey; filled triangles - GEMS, open squares - GOODS, 
stars - HUDF. The solid line indicates the best-fit power law. 



Number^ 


Survey 


yPHOT 




„PHOT c 
' e 






(AB mags) 


(") 


(") 


G75 


HUDF 


27.07 ±0.04 


0.09 


0.15 


G76 


HUDF 


28.21 ±0.10 


0.35 


0.10 


G78 


HUDF 


26.09 ±0.01 


0.19 


0.13 


G83 


HUDF 


26.23 ±0.01 


0.35 


0.14 


G90 


HUDF 


27.04 ±0.03 


0.50 


0.14 



^ Index 

^ Distance between ACS and ground-based centroids 

Half-light radius computed by PHOT (not reported for 
LAEs without SExtractor detections) 
* This table is only a stub. A manuscript 
with complete tables is available at 
http: / /www. nicholasbond.com/Bond0413.pdf 



z — 3.1 and, considering only z2EWcomplete, 108 LAEs 
at z = 2.1. When an object is covered by multiple sur- 
veys, we use the cutout from the deeper survey. We 
note that there are two objects in z2EWcomplete (G152 
and G180) ostensibly covered by the GOODS survey that 
have obvious defects in the Vgoe-band images. Both of 
these objects also have GEMS coverage, so we instead 
analyze their defect-free GEMS cutouts. 




(") 



Fig. 10. — Fractional uncertainty in the half-light radius as a 
function of half-light radius for a series of simulated Gaussian- 
shaped light profiles with constant continuum magnitudes. Each 
curve represents a different continuum magnitude with, from top 
to bottom, Veoe = 24.2, 25.0, 25.7, 26.5, and 27.2. Each cross 
symbol indicates the median over 1000 simulations. 

magnitudes because of the difficulties associated with 
measuring the light centroid (the center of the aperture) 
on a discrete grid. This effect does appear in our half- 
light radius simulations, but it results in only a ~ 4% 
overestimate of for marginally resolved sources with 
Vgoe =27. The effect is even smaller for objects that are 
more extended and/or brighter than these limits. 

4. LAE MORPHOLOGY EVOLUTION BETWEEN Z = 3.1 
AND Z = 2.1 

Between the GEMS, GOODS, and HUDF surveys, 
there is HST/ ACS coverage for a total of 171 LAEs at 



4.1. UV continuum photometric centroid 

The UV continuum emission in an LAE's host galaxy 
travels directly to us from the young stars, but Lya neb- 
ular emission can resonantly scatter to large distances 
from the original starburst, depending on the distribu- 
tion of neutral hydrogen surrounding the host galaxy. In 
this circumstance, we might expect to see an offset be- 
tween the emission-line centroid, measured in a narrow- 
band filter from the ground, and the rest-frame UV cen- 
troid, measured in the Vgoe HST/ ACS cutouts. 

Using the procedure described in Section I3.2[ we 
compute the Vgoe-band centroids of all of the objects 
with HST coverage. There were five objects, one in 
zSCiarduUo (C34) and four in z2EWcomplete (G168, 
G206, G218, and G222), for which a centroid could 
not be computed because SExtractor did not detect any 
sources within rgd- A visual inspection of their cutouts 
reveals no evidence for likely counterparts in C34, G206, 
or G222, but the remaining two have extended sources 
just outside rgei- For the computation of fixed-aperture 
magnitudes and half-light radii, we set the centroids of 
these objects to their ground-based, narrow-band posi- 
tion (the cutout center). 

In Figure 111! we plot the distribution of measured 
offsets between the Vgoe-band continuum centroids 
and ground-based emission-line positions in both the 
z2EWcomplete catalog and the combined z = 3.1 LAE 
catalog. We fit two-dimensional Gaussians to the distri- 
butions of centroid offsets and find a = 0.24 ± 0.04 and 
a — 0.20 ± 0.03 for z = 2.1 and z = 3.1, respectively. 
This amount of scatter is consistent with the expected 
0'.'2 — 0'.'3 astrometric uncertaintes in the ground-based 
narrow-band surveys, so we find no evidence for a physi- 
cal offset between the emission-line and continuum light 
distributions at the ^ 0'.'2 level (^ 1.5 kpc at z = 2 — 3) 
at either redshift. 
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TABLE 3 
z = 3.1 LAE Properties 



Number'' 


Survey 




■5^ 


yPHOT 




„PHOT d 
' e 










(AB mags) 


(") 


(") 


C15 


GOODS 


3:33:03.233 


-27:50:48.260 


24.99 ±0.05 


0.49 


0.10 


C22 


GOODS 


3:32:38.849 


-27:41:44.061 


24.80 ±0.03 


0.23 


0.20 


C34 


GOODS 


3:32:17.457 


-27:42:45.900 








C49 


GOODS 


3:32:33.117 


-27:54:19.552 


26.32 ±0.08 


0.28 


0.23 


C52 


GOODS 


3:32:54.834 


-27:46:40.132 


25.48 ±0.07 


0.14 


0.13 



Index from table 2 of CiarduUo et al. 2010 
^ Position of ACS centroid (set to ground-based position when there are no SExtractor 
deteetions) 

Distance between ACS and ground-based centroids 

Half-light radius computed by PHOT (not reported for LAEs without SExtractor de- 
tections) 

* This table is only a stub. A manuscript with complete tables is available at 
http: / /www. nicholasbond.com/Bond0413.pdf 
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Fig. 11. — Normalized distributions of distances between the 
Veoe-band centroids and the ground-based narrow-band positions 
of LAEs at 2 = 2.1 (solid histogram) and z = 3.1 (dashed). The fre- 
quency distribution for the best-fit two-dimensional Gaussians to 
each centroid offset distribution are also plotted as solid (cr = 0'.'24) 
and dashed (cr = 0'.'20) curves. The dispersions of these Guassians 
are consistent with the uncertainties in the ground-based astrom- 
etry (0'.'2 - 0'.'3). 

4.2. Fixed- aperture photometric properties 

The fixed-aperture half-light radius measures the size 
of the LAE "system;" that is, the size of the combined 
light distribution within rgoi of the rest-UV continuum 
centroid. If, for example, an LAE originates from a 
pair of merging galaxies with comparable UV continuum 
brightness, then this half-light radius would be approx- 
imately the distance between the galaxies. We give the 
fixed-aperture Vgoe magnitudes, half-light radii, and cen- 
troid offsets (relative to the catalog position) for the ob- 
jects in z2EWcomplete and zSCiarduUo in Tables [2] and 
[3l respectively. The corresponding tables for zSGronwall 
can be found in B09. 

In Figure fT2l we plot the distribution of fixed- aperture 
half-light radii for z2EWcomplete and the combined 
z ^ 3.1 LAE samples. Where a large fraction of the 
z — 3.1 LAEs are near the resolution limit, the major- 
ity z = 2.1 LAE systems are resolved, with a tail ex- 



tending to larger half-light radii. The median half-light 
radii are r™'^^ = 1.41 kpc and r™°'^ = 0.97 kpc for 
z2EWcomplete and the combined z = 3.1 sample, re- 
spectively. Using the Kolmogorov-Smirnov (K-S) test 
and making the null hypothesis that the two samples are 
drawn from the same half-light radius distribution, we 
find a P-value of 0.0003. If we exclude those objects 
in z2EWcomplete classified as DEx in a visual inspec- 
tion (and therefore less likely to be high-redshift LAEs, 
see Section [3?T|l . we find P = 0.004. Finally, even if we 
exclude all objects with large half-light radii, requiring 
^PHOT ^ 2.5 kpc, we find P = 0.05. 

We note that the Vgoe filter is centered on rest-frame 
wavelengths of ~ 1500 A and 2000 A at z = 3.1 and 
z = 2.1, respectively. Although emission at both wave- 
lengths is dominated by UV radiation from young stars, 
there may still be a difference in the apparent sizes of 
LAEs between these two parts of the spectrum. To 
test for this, we compared the sizes of 18 z = 2.1 
LAEs in the observed-frame B- and V-band imaging from 
the GOODS survey. We found the B-band (rest-frame 

1500 A) sizes to be 13% larger on average. This may 
be due to diffuse Lya emission leaking into the B-band 
filter, but it is difficult to tell from broadband imag- 
ing alone. Regardless, these differences only act to in- 
crease the size evolution we measure between z = 2.1 and 
z = 3.1. We therefore conclude that, for LAEs selected 
with the same rest-frame equivalent width and Lya lumi- 
nosity cutoffs, the systems at z = 2.1 are systematically 
larger than those at z — 3.1. 

In Figure [T51 we show the dependence of fixed-aperture 
half-light radius on UV continuum magnitude for both 
z2EWcomplete and the combined z = 3.1 samples. 
To make the comparison more direct, we have added 
0.77 mag to the z — 2.1 Veoe-band magnitudes, corre- 
sponding to the cosmological dimming that would occur 
if they were seen at z = 3.1. At both redshifts, there 
is a great deal of scatter that is largely independent of 
continuum brightness and much larger than the observa- 
tional uncertainties (which are cr^/r 0.1 at Vgoe = 26, 
see Figure ITOl) . In order to test for a correlation between 
UV continuum magnitude and half-light radius, we have 
divided the data into five bins in continuum magnitude, 
computed the mean half-light radius within each bin, and 
estimated the uncertainty on this mean using 5000 boot- 
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TABLE 4 
2 = 2.1 LAE Component Properties 



Number^ 


Component'^ 


Survey 








6*° 


„SE f 
' e 


SFR(UV) 










(") 


(AB mags) 


n 


(") 


(Me/yr) 


G75 


1 


HUDF 


27.41 ± 0.02 


0.03 


0.47 


-17.30 


0.14 


0.37 


G76 


1 


HUDF 


28.29 ± 0.04 


0.27 


0.61 


-11.40 


0.08 


0.16 


G78 


1 


HUDF 


26.49 ± 0.01 


0.17 


0.82 


51.80 


0.08 


0.85 




2 


HUDF 


27.42 ± 0.02 


0.21 


0.70 


42.20 


0.09 


0.36 


G83 


1 


HUDF 


26.50 ±0.01 


0.37 


0.93 


-46.60 


0.07 


0.85 



^ Index from table 2 of Gronwall et al. 2007 
^ Component number 

Distance from ground-based Lya position 

Isophotal axis ratio computed by SExtractor 

Isophotal position angle computed by SExtractor 
* Half-light radius computed by SExtractor 

This table is only a stub. A manuscript with complete tables is available at 

http: / / www.mcholasbond.com/Bond0413.pdf 



strap simulations. The resulting means and their uncer- 
tainties are plotted as large points in Figure [T3l 

There is no evidence for a correlation between r^^'^'^ 
and Vgoe at z — 2.1 - a constant value of = 1.49 kpc 
is a good fit at — 2.6. However, the best-fit con- 
stant value to the combined z = 3.1 LAE sample, 
j„PHOT _ gg j^p^^ jg ^ significantly poorer fit to the data 

(x^ ~ 11.5) than a model where size decreases logarith- 
mically with continuum flux (x^ = 2.15). The best-fit 
two-parameter model to the z = 3.1 sample is 

r^^^^ikpc) = -0.21V606 + 28.8. (3) 

It is worth noting that this relationship only deviates 
significantly from the z = 2.1 data at faint continuum 
magnitudes, Vgoe ^ 26.5, suggesting that this is where 
most of the size evolution has occurred between z = 3.1 
and z — 2.1. 

We show a similar plot in Figure [T31 which gives the 
relationship between fixed-aperture half-light radius and 
rest-frame equivalent width. We find no evidence for 
a correlation between equivalent width and half-light ra- 
dius at either redshift. This is in line with what is seen in 
a sample of Lyman break galaxies at 2.5 < z < 3.5, which 
exhibits little dependence of the UV morphol ogy on the 
prese nce or strength of Lya emission (jPentericci et al.l 

[2oTnh . 

4.3. Properties of photometric components 

In the terminology of this paper, a photometric "com- 
ponent" is any contiguous source within rgoi of the 
ground-based Lya centroid of an LAE, where sources 
are identifie d in the Veoe HST images using SExtractor 
(see Section IX^ . We give the brightness, ellipticity, po- 
sition angle and observed half-light radius (rf^) of each 
LAE component (as computed by SExtractor) for the 
z2EWcomplete and z3Ciardullo sample in Tables |3] and 
O respectively. We quote half-light radii as computed 
by SExtractor rather than in fixed apertures in order to 
avoid blending effects in multi-component systems. 

Photometric components can be galaxies themselves or 
individual star-forming regions within a larger galaxy - it 
is difficult to distinguish these two possibilities from mor- 
phology alone for very compact objects like LAEs. Fur- 
thermore, objects with multiple photometric components 
can be chance coincidences with low-redshift galaxies. In 
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Fig. 12. — Distributions of fixed-aperture, observed half-light 
radii for LAE systems at 2; = 2.1 (solid and dotted) and 2 = 
3.1 (dashed). The dotted histogram includes all objects in 
z2EWcomplete with HST coverage and the s olid histogram ex- 
cludes objects classified as DEx (see Section 13.11 1. The verti- 
cal lines indicate the approximate resolution limit of the V-band 
HST images at each redshift. We include only LAE systems with 
Veoe < 28. 

z2EWcomplete, excluding objects classified as DEx, 19 
of 108 objects in the sample have more than one pho- 
tometric component. This is consistent with the 23 of 
171 objects (13%) with multiple components at z = 3.1, 
but it is also consistent with our estimated maximum 
rate of interlopers (~ 14%) in the selection region of 
the z = 2.1 LAEs (based on the density of background 
sources in the cutouts, see Section [321) • If the majority 
of multi-component systems in our LAE samples are con- 
taminated by low-redshift interlopers, the distribution of 
component separations should be consistent with what 
we would see when comparing the LAE central compo- 
nents to a population with a random angular distribu- 
tion. To simulate this, we keep the brightest component 
in each multi-component system and randomly place an- 
other component within the selection circle. 

We show the normalized distribution of separations for 
the multi-component systems in our data and 10'* mock 
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TABLE 5 

z = 3.1 LAE Component Properties 



Number'' 


Component'' Survey 


a 


6 


ySE 






rfE t SFR(UV) 












(") (AB mags) 


n 


(") (Mo/yr) 



C15 


1 


GOODS 


3:33:03.232 


-27:50:48.267 


25.12 ±0.01 


0.44 


0.56 


-67.20 


0.10 


3.02 


C22 


1 


GOODS 


3:32:38.850 


-27:41:44.056 


25.44 ±0.01 


0.22 


0.32 


-50.80 


0.20 


2.24 


C49 


1 


GOODS 


3:32:33.108 


-27:54:19.636 


26.90 ±0.03 


0.21 


0.68 


-57.40 


0.10 


0.58 




2 


GOODS 


3:32:33.135 


-27:54:19.354 


27.73 ±0.07 


0.43 


0.83 


-21.90 


0.08 


0.27 


C52 


1 


GOODS 


3:32:54.836 


-27:46:40.134 


25.75 ±0.02 


0.11 


0.82 


84.30 


0.11 


1.68 



^ Index from table 2 of CiarduUo et al. 2010 
^ Component number 

Distanee from ground-based Lya position 

Isophotal axis ratio computed by SExtractor 
° Isophotal position angle computed by SExtractor 
' Half-light radius computed by SExtractor 

This table is only a stub. A manuscript with complete tables is available at http://www.nicholasbond.com/Bond0413.pdf 
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Fig. 13. — Fixed-aperture, rest-UV half-light radius plotted ver- 
sus rest-UV continuum magnitude in the full sample of LAEs with 
SExtractor detections, including objects at z = 3.1 (solid trian- 
gles) and 2 = 2.1 (open squares). We have added 0.77 mag to the 
2 = 2.1 Veoe-band magnitudes, corresponding to the cosmological 
dimming that would occur if they were seen at z = 3.1. Small 
points indicate individual LAEs, while large points with error bars 
indicate the median half-light radii in bins of AVeoe = 0.6. Un- 
certainties on the median are each computed from 5000 bootstrap 
simulations. The dotted line indicates the approximate resolution 
limit of the Veoe-band HST images. 

multi-component systems in Figure [T51 At z = 3.1, the 
real multi-component systems have a median separation 
of 0'.'41, while the same statistic is 0'.'58 for the mock 
systems. With the null hypothesis that the two sam- 
ples are drawn from the same parent distribution, the 
K-S test yields P = 0.029, suggesting that the major- 
ity of the components in the multi-component systems 
in zSCiardullo and zSGronwall are associated with LAEs 
at z = 3.1. By contrast, we cannot exclude the null hy- 
pothesis for z2EWcomplete (P > 0.1), despite a smaller 
median separation in the observed systems (O'.'Sl) than 
in the mock systems (0'.'58). Although this may mean 
that the multi-component systems in z2EWcomplete are 
heavily contaminated by interlopers, this result would 
also be consistent with an increase in the median sep- 
aration between components in multi-component LAEs 
from z = 3.1 to z = 2.1. 



Fig. 14. — Same as Figure (TS] except plotting fixed-aperture, 
rest-UV half-light radius as a function of rest-frame Lya equiva- 
lent width, where EW (Ly a) measurements are t aken from GulO, 
ICiarduUo eTall | [20Tll) . and lGronwall et al.lll200W) . The dotted line 
indicates the approximate resolution limit and the dashed line in- 
dicates the approximate rest-frame EW limit of the LAE surveys. 

Another way to discern the nature of the multi- 
component systems is to compare their component size 
distribution to that of single-component systems. If low- 
redshift interlopers are common in the multi-component 
systems, these components might also have larger half- 
light radii than the average LAE. We show the compo- 
nent size distributions in single- and multi-component 
systems in Figure [1^1 The K-S test reveals no evidence 
that the single- and multi-component systems are drawn 
from a different parent distribution at either z — 2.1 or 
z — 3.1 {P > 0.1). It is therefore unlikely that the multi- 
component systems at either redshift are dominated by 
low-redshift interlopers, unless their size distribution is 
very similar to that of high-redshift LAEs. 

In B09, we showed that the half-light radius mea- 
sured by SExtractor was systematically underestimated 
for components detected at S/N< 30. For the shallowest 
survey used here (GEMS), this corresponds to Vgoe — 
26.2, so we construct subsamples of components brighter 
than this limit (including those in multi-component sys- 
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Fig. 15. — Normalized distributions of separations between com- 
ponents in multi-component LAE systems. When there are more 
than two components, we count as one every pair that includes the 
brightest component. The solid and dotted histograms are the dis- 
tributions for the 2 = 2.1 and z = 3.1 LAE samples, respectively. 
The dot-dashed and dashed histograms give the distributions we 
would get (for z = 2.1 and z = 3.1, respectively) if the fainter com- 
ponent of each pair were randomly distributed within the selection 
circle (tscI = 0'.'65 at 2 = 2.1 and 0'.'6 at 2 = 3.1). 



z = 3.1 




r^^ (kpc) 




(kpc) 



Fig. 16. — Normalized SExtractor half-light radius distribution 
for components in single-component systems (solid) and multi- 
component systems (dotted) for the z = 3.1 (top panel) and 
2 = 2.1 (bottom) LAE samples. We only plot LAE components 
with V606 < 26.2. 

terns). At z — 3.1, the remaining components have 
med(r'f ^) = 0.79 kpc, compared to med(rf — 0.95 kpc 
at z = 2.1. With the null hypothesis that the two sam- 
ples are drawn from the same parent distribution, the 
K-S test gives P = 0.07, again consistent with an in- 
crease in the size of the typical LAE from z = 2.1 to 
z = 3.1. 

4.4. Median morphological properties of LAE 
subsamples 

iGuaita et al.l ()2011[ ) selected a series of subsamples 
from the z2Guaita LAEs, separately dividing the sam- 
ple by Lya equivalent width, rest-frame UV luminos- 
ity, 3.6 /zm brightness, and B— R color. We present the 



morphological properties of these subsamples in Table [51 
In addition to median half-light radius, we give the in- 
terquartile range (Q3-Q1) of the half-light radii. Since 
measurement errors are subdominant to intrinsic scatter 
in the LAE sizes (see Section 14. 2|) , the latter quantity 
will give an indication of the morphological heterogene- 
ity within the subsample. Uncertainties on each quantity 
are estimated from 10^ bootstrap simulations. 

For comparison, we also create a set of subsamples at 
z = 3.1. We use the same Lya equivalent width selection 
criterion as at z = 2.1, but correct the rest- frame UV 
luminosity and 3.6 /im brightness selection criteria for 
the difference in luminosity distance. Finally, for the 
rest-UV color, we use V— R instead of B— R, attempting 
to match the rest-frame wavelengths of the bluer band. 
The properties of each subsample are shown in Table [T] 

We find no significant size difference between the high- 
and low-equivalent width subsamples at either redshift, 
consistent with the object-by-object results shown in Fig- 
ure [141 The UV-bright subsamples are larger in both size 
and size spread than the UV-faint subsamples at both 
redshifts, but the difference is greater and more statis- 
tically significant at z = 2.1, with Ar^ — 0.36 ± 0.13. 
A statistically significant difference only appears when 
we include all of z2Guaita - when the subsamples are re- 
stricted to objects in z2EWcomplete, the half-light radius 
difference between UV-bright and UV-faint subsamples 
decreases to Are = 0.17 ± 0.18. 

The other two pairs of subsamples, those selected by 
3.6 fj,m flux and B— R color, also show a difference be- 
tween their median half-light radii at z — 2.1. In general, 
the subsample with the larger median size has properties 
consistent with what we would expect from objects with 
a dusty, older stellar population underlying the recent 
star formation. The brightness at 3.6 /im, in particular, 
is an indicator of the total stellar mass of an LAE and will 
be large when a galaxy has undergone previous bursts of 
star formation. At z = 3.1, we still see a significant size 
difference between the red and blue LAEs, but the frac- 
tion of red LAEs is smaller (11% as compared to 17% at 
z = 2.1), indicating that the dusty objects, although still 
larger in size at z = 3.1, make up a smaller fraction of the 
LAE population. We see no significant difference in size 
between the IRAC-faint and IRAC-bright subsamples at 
z = 3.1, and the fraction of IRAC-bright objects is again 
smaller (26%, as compared to 44%). 

All of the subsamples exhibit a size spread, as indicated 
by the interquartile range of the half-light radius distri- 
bution, that is considerably larger than that expected 
from the measurement uncertainties alone. The small- 
est spread (0.5 kpc) is seen in the UV-faint subsample, 
but its fractional spread (IQR[re]/med[re]) is compara- 
ble to that of the other subsamples. The most noticeable 
outlier is the IRAC-bright subsample, which exhibits a 
size spread of 1.2 kpc, suggesting an unusual amount of 
heterogeneity in this subsample. A comparable spread is 
not seen in the z = 3.1 IRAC-bright subsample. 

5. DISCUSSION 

The case for evolution in the LAE popul ation between 
z ^ 3 and z ~ 2 is becoming quite strong. iNilsson et al.l 
(|2009D . comparing samples of LAEs at z ~ 3 and z = 
2.25, claim an increase in the AGN fraction and the UV- 
to-Lya SFR ratio, as well as a narrower EW distribution 



10 



Bond et al. 



TABLE 6 
Sizes of z=2.1 LAE Subsamples 



Subsample^ 


Selection 


Number of LAEs'' 


Median Size"^ 


Size Spread'* 








(kpc) 


(kpc) 


UV-faint 


R > 25.5 


55 


1.13 ±0.09 


0.53 ±0.16 


UV-bright 


R < 25.5 


78 


1.49 ± 0.09 


0.84 ±0.13 


IRAC-faint 


h.e < 0.57 M Jy 


38 


1.10 ±0.11 


0.66 ±0.13 


IRAC-bright 


/3.6 > 0.57 M Jy 


30 


1.48 ± 0.27 


1.24 ±0.18 


low-EW 


EW < 66 A 


39 


1.23 ± 0.14 


0.69 ±0.12 


high-EW 


EW > 66 A 


31 


1.12 ±0.10 


0.64 ±0.21 


red-LAE 


B -R>0.5 


19 


1.64 ±0.21 


0.83 ±0.23 


blue-LAE 


B~ R<0.5 


114 


1.27 ±0.10 


0.77 ±0.09 



LAE subsamples defined in Guaita et al. 2010 
^ Includes only LAEs with HST/ ACS coverage 

Median fixed-aperture half-light radius within 0''65 
"* Interquartile range of fixed-aperture half-light radius 



at z = 2.25. Analyzing the sample studied in this pa- 
per, GulO confirmed an increase in UV-to-Lya SFR to- 
ward z ~ 2, but found no evidence for an increase in the 
AGN fraction aX z = 2.1. Also using the GulO sample, 
ICiardullo et al.l (j201lD find a ~ 50% decrease in the num- 
ber density for sources with L > 1.5 x 10''^ erg s~^, as 
well as a decrease in both the scale length of the equiva- 
lent width distribution and the characteristic luminosity 
(L*) of the luminosity function at z — 2.1. 

There is also evidence that the LAE po pulation is be- 
coming more heterogeneous with time. iNilsson et al.l 
(|2)11) found increased variety in the SEDs of LAEs at 
z = 2.25 as compared to higher-redshift samples, with 
only 15% of the galaxies being consistent with a single, 
young population of stars. Furthermore, GulO find a bi- 
modality in the B— R colors of the UV-bright {R < 25) 
objects in their sample at z = 2.1, possibly indicating 
the presence of a sub-population of LAEs for which a 
significant quantity of dust is present. The evidence for 
heterogen eity in the z2Guaita sample was further demon- 
strated in lGuaita et al.l ()2011f ). where they looked at the 
SED properties of a series of subsamples s elec ted using 
various photometric cuts (see also. Section l44l) . In Fig- 
ures [T71[TH1 andHni we show the relationship between the 
derived SED properties and the median half-light radii 
of the z2Guaita subsamplefl In all cases, we see a wide 
range of SED properties, as well as a clear correlation be- 
tween the LAE size and the properties derived from its 
SED. LAEs are found to be larger for galaxies with higher 
stellar mass, higher dust obscuration, an d higher star 
formation rate (averaged over 100 Myr, see lGuaita et al.l 
120 lit ). These results a re broadly consistent with th e nu- 
merical simulations of IShimizu fc Umemural ()20100 . who 
predict that size, stellar mass, and star formation rate 
will correlate positively with the mass-weighted age of 
LAEs. 

The apparent heterogeneity of LAEs at z^ 2 might 
indicate the presence of a sub-population of LAEs that 
are more massive and more evolved. Such objects would 
normally have their Lya emission extinguished by dust, 
but if the galaxies are accreting a relatively pristine 
subhalo, the Lya emission may be originating from a 
dust-free star-forming region sufficiently separated from 

^ Note that the plotted points are not all independent, as there 
is overlap in the subsamples 



the parent halo that the line emission could escape 
(HE imizu fc Umemural I2010D . Alternatively, the dust 
in these objects may be concentrated in dense clumps, 
with the Lya photons resonantly scattering o ff of the 
clunip surface and escaping the galaxy (e.g., iNeufeldl 
[iMl|Hansen & Oh"2006'; 'Finkelstein et ay|2009D. In the 
simulations of [Shimizu & Umcmura (20ic|, only ~ 30% 
of LAEs at z = 3.1 are evolved galaxies experiencing de- 
layed accretion of a subhalo onto a parent halo rather 
than galaxies undergoing their first major burst of star 
formation. By z — 2.1, however, their models show the 
majority of LAEs 70%) are in the former category. 
This is qualitatively consistent with the increase in size 
and the broader range of morphological properties that 
we see at z = 2.1, as well as the increase in dust redden- 
ing seen bv lGuaita et"all (poTTI ). 

More work needs to be done before we can properly 
quantify the size evolution of LAEs through cosmic time. 
Optimally, we would like to measure size evolution at 
z < 1.5 and z > 4 in order to determine whether the 
evolution seen here constitutes a sudden change in the 
LAE population or a more gradual evolution with red- 
shift. iTaniguchi et al.l (|2009| 1 studied the sizes of a sam- 
ple of z = 5.7 LAEs, fitting a PSF-convolved model 
to a stack of 43 LAEs and found a best-fit half-light 
radius of 0.76 kpc. This can be roughly compared to 
the median GALFIT-derived size of z = 3.1 LA Es in 
zSGronwall, rf = 70 kpc (IGronwall et al.] [2010l). but 
with the caveat that ITaniguchi et al.l ()2009f ) were work- 
ing with a sample that was considerably brighter in Lya 
(Li^va > 7 X 1 0-^^ erg s~^ iMuravama et al.l [20071 ) and 
IGronwaU et all (12010 ') were restricting themselves to in- 
dividual LAEs with S/N> 30 in the HST/ ACS Vgoe im- 
ages. 

It is possible that LAEs will eventually be found to re- 
fiect the approximately H~^(z) size evolution of the over- 
all galaxy population (' Ferguson et al.l [20041 ) . but there 
are good reasons to think that this will not be the case. 
LAEs can appear visible to the observer only so long as 
the Lya radiation is able to escape the galaxy without be- 
ing absorbed by dust. As such, they likely exist for only 
a short time after galaxy-scale star formation has turned 
on and before enough dust is formed to extinguish the 
Lya. Simulations suggest that this time p eriod could be 
shorter than 3 x 10^ years (e.g.. lMori fc Um emura 2006^ . 
If so, then LAEs are only a single snapshot in the history 
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TABLE 7 
Sizes of z=3.1 LAE Subsamples 



Subsample 


Selection 


Number of LAEs^ 


Median Size"^ 


Size Spread'^ 








(kpc) 


(kpc) 


UV-faint 


R > 26.3 


80 


0.97 ± 0.06 


0.55 ±0.07 


UV-bright 


R < 26.3 


56 


1.09 ±0.09 


0.75 ±0.12 


IRAC-faint 


fs.a < 0.3 At Jy 


31 


1.24 ±0.14 


0.65 ±0.08 


IRAC-bright 


/a. 6 > 0.3 fi Jy 


11 


1.24 ±0.24 


0.73 ±0.23 


low-EW 


EW < 66 A 


85 


1.06 ± 0.08 


0.67 ±0.10 


high-EW 


EW > 66 A 


51 


0.98 ± 0.06 


0.59 ±0.11 


red-LAE 


y - i? > 0.3 


13 


1.67 ±0.36 


1.40 ±0.30 


blue-LAE 


V - R<0.3 


123 


1.00 ±0.05 


0.58 ± 0.07 



" Includes only LAEs with HST/ ACS coverage 
Median fixed-aperture half-light radius within 0'.'6 
Interquartile range of fixed-aperture half-light radius 



UV-bright 

UV-fainl 

IRAC-bright 

IRAC-faint 

recJ-LAE 

blue-LAE 

low-EW 

high-EW 



UV-bhght 

UV-faint 

IRAC-bhght 

IRAC-faint 

recJ-LAE 

blue-LAE| 

low-EW 
high-EW 



-<) 1 



log Stellar Mass (M^^j^) 



E(B-V) 



Fig. 17. — Median fixed-aperture half-light radius versus stellar 
mass in subsamples of 2 = 2.1 LAEs. Stellar masses are taken from 
an SED median stacking analysis in Guaita et al. 2010. 



Fig. 18. — Same as Figure [TTl but plotting fixed-aperture half- 
light radius as a function of dust reddening, E(B — V). 



of a galaxy's formation and their size evolution should be 
much less steep than that of the overall galaxy popula- 
tion. If, on the other hand, Lya emission is able to escape 
from a large fraction of galaxies with previous genera- 
tions of stars alrea dy in place, they should more closely 
trace the Ferguso n et al.l (|2004[ ) law. The heterogeneity 
seen in the present samples suggests that the LAE pop- 
ulation contains both types of objects, but more work 
is needed to elucidate whether this division is an actual 
bimodality or simply a continuous range of properties. 
The z — 2.1 LAE sample, in particular, may contain up 
to ~ 15% contamination from low-redshift galaxies (see 
Section 13.11) - spectroscopic follow-up is needed to ac- 
curately estimate the contamination fraction and isolate 
the types of objects that contribute to the contamina- 
tion. Furthermore, analysis of the deep rest-frame opti- 
cal (observed- NIR) imaging obtained as part of th e Wide 
Field Camera 3 Early Release Science (Windhor st et al.l 
I2010D and the Cosmic Assembly Near-infrared Deep Ex- 
tragalactic Legacy Survey would allow us to determine 
where most of the stellar mass in these objects lies and 
whether it coincides spatially with the rest-UV emission 
from the young stars. 
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Fig. 19. — Same as Figure (TT] but plotting fixed-aperture half- 
light radius as a function of star formation rate. 
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Fig. 1. — Cutouts of 2 = 2.1 LAEs extracted from the GEMS survey images. We mark components (SExtractor detections within 0'.'65 
of the cutout center) with arrows. Numbers underneath the panels are the corresponding LAE indices from GulO. 



14 



Bond et al. 









■' 




iT r --w-r- : 

■ ■ - ' - ■ 




15 




47 
































































^ .; 




. . ' 
































. ■■ ;, A -ii; 

■ - . v-rv. 
















6 


17 


39 














t - 














































I B 40 H pa 

2 




P3 



1. 





















ei 



St 

> J.' 



SB 



■ i^-v^ 'V 



S4 





t 











B3 



ill: 



64 



63 



r r |- 



J L 



m 



Fig. 1. — (cont.) 
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Fig. 3.— 



Same as Figure [T] but for LAEs covered by the GOODS Vgoe-band images. 
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Same as Figure [2] but for LAEs covered by the GOODS Vgoe-band images. 
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Same as Figure [T] but for LAEs covered by the HUDF Vsoe-band images. 
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